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In this Letter we study a new class of synthetic materials in which the refractive index satisfies 
a special symmetry, n(—x) = — n*{x), which we term antisymmetric parity-time (AVT) systems. 
Unlike PT-symmetric systems which require balanced gain and loss, i.e. n(— x) — n*(x), AVT 
systems consist of balanced positive and negative index materials (NIMs). Despite the seemingly 
PT-symmetric optical potential ^(a;) = n(x) 2 uj 2 / c 2 , such systems are not invariant under combined 
VT operation due to the discontinuity of the spatial derivative of the wavefunction. We show that 
AVT systems display intriguing properties such as spontaneous phase transition of the scattering 
matrix, bidirectional invisibility, and a continuous lasing spectrum. 
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Parity-time {VT) symmetric systems have attracted 
considerable interests since the original work by Ben- 
der and colleagues [TH3] , both in the fields of quan- 
tum mechanics and photonics. In 'PT-symmetric quan- 
tum mechanics the Schrodinger equation does not have 
time-reversal symmetry but is invariant under combined 
parity-time operation. The corresponding Hamiltonian 
is non-Hermitian, yet it displays regions of real energy 
spectrum and was speculated as an alternative to canon- 
ical quantum mechanics. The idea was introduced into 
photonics by Christodoulides and coworkers |H [5], uti- 
lizing balanced gain and loss (n(— x) = n*(x)) to satisfy 
the PT-symmetry. Many intriguing optical phenomena 
have since been predicted and observed, such as double 
refraction [5], power oscillations [MI], coexistence of co- 
herent perfect absorption [9] and lasing states [10l [11] , 
spontaneous symmetry breaking of the scattering matrix 
[llj , and unidirectional transmission resonances [121 113) . 

Another field that has attracted tremendous research 
interests is negative index materials (NIM) [HHT7] . In 
such materials the refractive index is negative over some 
frequency range, achieved by using engineered resonances 
of nanostructures in metamaterials. NIMs have been pro- 
posed to realize subwavelength imaging (i.e. "superlense" 
[13]) and dramatically change light-matter interaction 
via altered local density of optical states [TH] . NIMs are 
strongly absorptive, but recent advances with imbedded 
active media have demonstrated reduced intrinsic loss 
and even net gain |19j . 

In this Letter we study a class of new systems bridg- 
ing NIMs and PT-symmetric photonics. Their refractive 
index is antisymmetric under combined VT operation, 
n(— x) = —n*(x), i.e. with balanced positive index ma- 
terials (PIMs) and NIMs; the imaginary part of n(x) is 
symmetric, which can be positive (loss), negative (gain), 
zero, or any complicated spatial function. Outside the 
system we assume that the refractive index is uniform, 
symmetric, and positive. We term such synthetic sys- 
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FIG. 1. (Color online) Schematic of a ID AVT photonic 
heterostructure, consisting of 2N layers that satisfy n(x) = 
—n(—x)*. 



terns antisymmetric parity-time (.4PT) systems, which 
display intriguing features such as bidirectional invisibil- 
ity, spontaneous phase transition of the S-matrix, and a 
continuous lasing spectrum. 

We base our discussion on the scalar wave equation for 
the electric field 



[V 2 + n{x) 2 (w 2 /c 2 )] E(x; uj) = 0, 



(1) 



which describes steady-state solutions in one-dimensional 
(ID) and two-dimensional (2D) systems. Henceforth we 
set c = 1 and assume that the corresponding /inim = 
— ImpimI if ^nim = — n piM- By nrs t gl anc e one may mis- 
judge that all the intriguing phenomena found in con- 
ventional PT-symmetric systems would survive since the 
optical potential V(x) = n(x) 2 ui 2 is still invariant under 
VT operation. We note that this is not true since the 
boundary condition at PIM and NIM interfaces is now 
different. Take a ID AVT heterostructures for example 
(see Fig. [I]), the electric field itself is still continuous at 
PIM and NIM interfaces, but its spatial derivation now 
satisfies gUl [H] 
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which changes abruptly due to the sign difference of /ipim 
and /xnim- Below We first analyze wave propagation and 
lasing in such ID AVT heterostructures, followed by a 
discussion of pseudo-yWT symmetry for wave propaga- 
tion in 2D with the paraxial approximation. 

In Ref. [11] the phase transition of the scattering ma- 
trix (S-matrix) in a PT-symmetric system is predicted 
based on the invariance of the system under combined 
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PT operations; We would not expect a similar phase 
transition in APT systems, since flipping the sign of the 
real part of the refractive index is not related to any sym- 
metry operation of the physical state. However, there do 
exist some special properties of the transmission coeffi- 
cient t and the left and right reflection coefficients tl , tr 
in a ID APT hctcrostructure: 



r L = r* R , Im[t] = 0. 



(3) 



To understand these properties, we start by noting one 
observation: By changing the refractive index of each 
layer in an arbitrary photonic hetereostructure to its neg- 
ative complex conjugate and flipping the sign the mag- 
netic pearmeability, i.e. u, — > —n*, \ii — > —fii, the trans- 
fer matrix M [22], defined by 



= M 



(4) 



becomes its complex conjugate at the same real fre- 
quency: 



(5) 



The field amplitudes A, B, C, and D are defined in Fig.[T] 
or more specifically, 

f Ae -in w(x+L/2) + Be in u,(x+L/2) ^ x < -L/2, 
h(x;u>) — < Q e -i no u(x-L/2) _|_ £) e in u(x-L/2) ^ x > L/2, 

(6) 

where no is the refractive index outside the heterostruc- 
ture and we have assumed fio = 1. The proof of |5]) is 
straightforward from the analytical expression of M: 



(7) 



obtained from the continuity of E(x; oj) and Eq. |2|. The 
matrices Dq and rrii are given by 



Dn = 



1 1 

n -n 



cos(niwAj) i^ 1 sin(rija;Aj) 
i— sin(n i wAj) cos(rija;A i ) 



(8) 



(9) 



where A$ is the width of the ith layer. Under the trans- 
formation rii — > — n* , /ij — > — /ij (i = 1, . . . , N), one finds 
TOi(u) — > m*(w) at a rea/ frequency and so does M(u>). 
Since M(u) determines the wave propagation, all related 
quantities, such as r^, ru, t or equivalently the S-matrix, 
become their complex conjugate under this transforma- 
tion. 

We refer to the left half of an APT system U and the 
right half V. Eq. (JsJ implies that the transfer matrices 
of U and V arc related by 



Mi 



•H=a[M [ 7 1 H]V, a= JV 



(10) 



If we define Mu((S) = ( 



"in mi2 
m2i m 2 2 



and the transfer matrix 



of the whole APT system M avt (uj) = ("i 11 ™P " 
then find through M A -p-j-(uj) — Mjj (u>) My (oj) that 



m u = |mii| - 
m' 22 = |m 22 | 2 - 
m' 12 = mi 2 TO 22 

The S-matrix defined by 
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(11) 
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can be easily expressed in terms of Mj^-pq-: 
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(13) 



from which we immediately find ([3]) using (11 1. 

The phase transition of the S-matrix can be inferred 
from the relations ^ , which suggest the parametrization 
of the S-matrix by three independent real quantities: t, 
a = Re[rt], and b = Im[ri], i.e. 



S = 



a + ib t 
t a — ib 



(14) 



which is pseudo-Hermitian |23j . i.e. = r/Sr) 1 with 
V = ( l o ) ■ The eigenvalues of the S-matrix are given by 



s ± =a± \/t 2 -b 2 , 



(15) 



which have two phases and the phase transitions occur 
at t = ±b. The scattering eigenstates ip±(u>) = yc± 
display a transition simultaneously: 



p± = ^ = l[ib±Vt 2 -b 2 ]. 
u± t 



(16) 



In one phase (\t\ > \b\) the intensities of the two in- 
cident beams in a scattering eigenstate are the same, 
i.e. \p±\ = 1 (see Fig. [2]). Thus we refer to this phase 
as the symmetric scattering phase and the other as the 
symmetry-broken phase (\t\ < \b\). In the symmetric 
phase s± are real, meaning that the symmetric inputs are 
either amplified or damped equally with no phase added 
during the scattering process. In the symmetric-broken 
phase the two scattering eigenstates have the same scat- 
tering strength \s\ (s + = si) but with \p + \ = \p_\~ 1 . 

Tuning to the phase transition points requires adjust- 
ing either the gain/loss parameter of the system, the 
frequency of incident beam, or scaling the system size. 
Given that it is challenging to maintain APT (or PT) 
symmetry in the first two approaches due to material dis- 
persion, the third approach is probably most practical, 
i.e. by fabricating multiple scaled heterostructures and 
fixing the frequency of incident beams at a value that 
achieves the APT symmetry. 
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FIG. 2. (Color online) Phase transition of the S-matrix in 
an APT heterostructure. Top: An 8-layers example with 
Re [m] = 1.3, -2, -1.7, -3, Im [n(x)] = 0.04 (loss), and width 
At = 1.2, 0.996, 0.165, 0.531/mi (i = 1,2,3,4) for x > 0. 
Thick solid black lines in the middle and bottom panels show 
the transitions of the asymmetry factor \p\ and the scattering 
strength |s|. Thin solid green line in the middle panel indi- 
cates \t 2 — Im[rL] 2 | -1 , which approaches infinite at the phase 
transition points marked by the dashed vertical line. 



The phase transition discussed above is a general prop- 
erties of all ID APT systems, independent of whether 
the system has net gain or loss. By flipping the sign 
of Im[n(x)], i.e. changing local gain into loss and vise 
versa, the system merely undergoes a time reversal and 
the phase transitions happen at exactly the same loca- 
tions. 

There is only one exception which occurs when the 
local gain/loss is zero [24], i.e. Im [n(x)] = 0. In this 
case an APT heterostructure is always in the symmet- 
ric phase. More strikingly, the relations (|3| now take a 
special form: 



t = 1, r L = r R = 0, 



(17) 



i.e. an APT system becomes invisible, and it is indepen- 
dent of the complexity and size of the heterostructure and 
at what frequency the APT symmetry occurs. This phe- 
nomenon is robust upon a slight breakdown of the APT- 
symmetry or in the presence of a small Im [n(x)] =/= 
(see Fig. 3ja)). In addition, the invisibility is indepen- 
dent of whether the APT structure is standalone or inte- 
grated in a photonic environment, as long as the nearest 
elements have the same refractive index (see Fig. [3]^b)). 
If the APT symmetry can be maintained over a finite 
frequency range, a pulse transmitted within this fre- 
quency window will be exactly the same as the initial 
pulse, with no pulse distortion or shrinking/expansion. 
This phenomenon is independent of the propagation di- 
rection, in contrast to the one-way invisibility found in 
■pT-symmetric heterostructures (12) . 

The relations (JTYJ) can be treated as a generalization of 
the phenomenon of vanished reflection that happens at 




FIG. 3. (Color online) Top: Transmittance and reflectance 
in the APT photonic heterostructure shown in Fig. [2] but 
with Im [n(i)] = 10 -4 (loss). Thick black line, thin blue line, 
and dotted red line indicate T — \t\ 2 , R — \tl\ 2 ~ \rit.\ 2 , 
and Arg[i], respectively. They are very close to the values 
given by |17l, i.e. 1, 0, 0, when Im [n(x)] = 0. Bottom: 
Same as the top panel but Re [714] changes from -3 to 3. The 
oscillations in R and T have a single period determined by 
the combined width of the leftmost and rightmost layers, i.e. 
AuiL = 7rL/(2A4Re [714]) = 5.64, as if the central 6-layer 
APT structure is absent. 



the interface of two impedance matched PIM and NIM 
materials (see Ref. [TJ] for example). There is at least 
one such interface in an APT system, i.e. at x = 0, 
but multiple reflections occur at other interfaces between 
two NIMs, two PIMs, and a NIM and a PIM of different 



\n\. One way to prove (17) is from the transfer matrix 
M_apj-(uj) directly: 



M Avr {u) = D^ 1 [iPl^m,] D , i ± 0. 



Note that 



rri—i(u})mi(u>) = 1 



when rii is real. Therefore, 



m avt (uj) = Do 1 [n 



i=-N,. 



-2,2, ...,N mi] D = 



which implies (17) 



(18) 



(19) 



= 1, 
(20) 



Net gain has been demonstrated in NIMs by imbed- 
ding an active medium [19] . which opens the possibility 
of achieving lasing in metamaterials. Lasing modes are 
given by the poles of the S-matrix on the real frequency 
axis, which in general lead to discrete set of solutions con- 
taining the lasing frequency and the corresponding 
threshold = — lm[n(x)] > assuming a spatially 

uniform gain. Each lasing mode has its distinct inten- 
sity profile and roughly speaking the mode order m in- 
dicates the number of peaks inside the cavity. To deter- 
mine {ojl,t} for each mode, one can, for example, solve 
the two equations given by the real and imaginary part 



of m'22 in ID heterostructures (see Eq. (13)). However, 
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FIG. 4. (Color online) (a) Threshold value t(lul) for the 
continuous lasing spectrum in a 2-layer AVT heterostruc- 
ture. Each layer is 500nm thick and the refractive indices are 
±2 — it(ojl) at threshold. Squares indicate the discrete lasing 
solutions u)L = 1.735, 3.245, 4.786 in a uniform PIM cavity of 
the same length. Inset: Intensity profiles at uiL = 1.5, 3, 4.5 
marked by the red dots in the main figure. Shadowed areas 
indicate the cavity, (b) Number of lasing modes in the region 
fcl/o < 30 versus the width of one layer in (a). 



Eq. ( 11 ) shows that m! 2 2 is always real for an AVT sys- 
tem, which implies that there exist continuous regions of 
lol in which a w^-dependent r can be found. In other 
words, a lasing mode can appear anywhere in these fre- 
quency regions and a mode order m cannot be assigned. 
A 2-layer AVT structure in the low frequency regime is 
shown as an example in Fig. [4] We observe a reduced 
threshold as the lasing frequency increases and more os- 
cillations gradually appear in the intensity profile. If we 
consider the lasing modes in the corresponding PIM cav- 
ity of the same length and \n\, we find that its discrete 
lasing solutions lie exactly on the continuous threshold 
curve of the AVT structure (see Fig. [1] and Appendix Ah. 
This comparison also shows that coherent feedback does 
occur in a .4"PT-system based laser, even though the its 
spectrum is continuous. 

Unlike the bidirectional invisibility, the continuous las- 
ing spectrum is singular and breaks down if the AVT- 
symmetry is broken, which can be utilized as a sensitive 
measure of the quantities of interest that lead to the lat- 
ter. Here we consider one scenario where the the AVT- 
symmetry is broken due to a slight length mismatch. 
Consider a 2-layer cavity of length L similar to that stud- 
ied in Fig. [4] but with the PIM layer 5 longer than the 
NIM layer. We found that all lasing modes disappear 
except the ones at kL « mnL/(\n\S) (m = 1,2,...), 
whose thresholds are about the same when the AVT- 
symmetry is satisfied. Take L = 250/ito, S — 250 nm, 
and |n| = 2 for example, the corresponding wavelengths 
are A = m /im, which are are well separated and easy to 
achieve single-mode lasing. These modes originate from 
the resonances of the tiny section of length 5, which acts 
as an external cavity for frequency selection. In this ex- 
ample, the variation of S (or L) is enhanced by four times 
in the wavelength of the fundamental mode (m = 1) since 
A A = 4A<5/m, which can be easily measured. As a com- 
parison, the sensitivity to detect S is reduced by a factor 



of 6/L — 10~ 3 using a lasing mode of roughly the same 
wavelength in a uniform PIM cavity of length L, which 
also has a much denser spectrum to analyze. We note 
that the laser linewidths are comparable in the two sys- 
tems, since the thresholds of the corresponding modes 
are about the same as discussed. 

So far we have discussed the AVT symmetry with bal- 
anced PIMs and NIMs. One may attempt to realize a 
"pseudo-„4'PT" symmetry using only PIMs (or NIMs), 
satisfying n{x) — —n(~x)* but with /i(x) = fi(—x). 
One example in which such symmetry can be realized 
is wave propagation in 2D paraxial geometry 4, 5., with 
transverse index variation n(x) = tiq + Sn(x) satisfying 
|(5n(a;)| uq and Sn(x) = —5n(—x)*. The Helmholtz 
equation ([TJ in this case becomes [5] 



.d(j> 
dz 



dx 2 



+ 5n(x)ko 



= 0, 



(21) 



where 4>{x, z) is the slowly varying part of the electri- 
cal field E(x, z) = (f)(x, z)e lk ° z and x — y/2n k x is the 
scaled coordinate in the transverse direction. The trans- 
verse optical potential now is proportional to 8n(x) in- 
stead of index squared in the original Hcltmoltz equa- 
tion ([I]), and the intriguing phenomena discussed above 
disappear. The only exception happens when the sys- 
tem becomes equivalent to a conventional "PT-symmetric 
structure. The latter occurs, for example, if 5n(x) = 
Asm x+iB cos x (A, JJeR); shifting i by 7r/2 transforms 
6n(x) to Acosx — i-Bsini, satisfying n(x) — n(—x)*. 

In summary, we propose a new class of synthetic mate- 
rials which are antisymmetric under a combined parity- 
time operation, i.e. n(x) = —n(—x)*. AVT systems 
demonstrate interesting features such as bidirectional in- 
visibility, spontaneous phase transition of the S-matrix, 
and a continuous lasing spectrum. Properties of AVT 
systems in higher dimensions are under investigation and 
will be reported elsewhere. 
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Appendix A: Threshold correspondence in a simple 
AVT heterostructure and a PIM cavity 

The threshold t(ujl) of the continuous lasing spectrum 
in a simple 2-layer cavity discussed in the main text is 
given by the solution of the real equation 



cos a 



Im 



sin a\ 



IT + 



sina(cosa)* 



(Al) 



in which n r is the real part of the refractive index in the 
PIM, a = (n r — ir)ujLL/2, and L is the cavity length. 
In comparison, the threshold and the discrete lasing fre- 
quency in an uniform PIM cavity of the same length 
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are simultaneously determined by the following complex 
equation 



cos(2a) = i 
It implies that 



n r — it + 



1 



n r — it 



sm a cos a 



tana = —i(n r — it), — - 



(A2) 



(A3) 



or 



Re 



tana 
(tana)* 



By taking the real part of both sides of ( A2 1 after multi- 
plying them by (cos a)* /cos a, we find 



|cosa| + Rc 
= -Im 



tana 



(tan a)* 



sm a 



n r — it H I sina(cosa)* 

n r — it , 



, (A5) 



(A4) from which we recover (Al) using (A4) 
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